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ABSTRACT 

Context. In a previous work, three bright mid-infrared/radio sources were discovered in the nuclear region of startburst/AGN galaxy 
NGC1365. 

Aims. The present study aims at confirming that these sources are indeed young and massive "embedded" clusters, and at deriving 
their physical parameters, such as extinction, age and mass. 

Methods. Using ISAAC and VISIR at the VLT we obtained maps and low resolution spectra in the near- and mid-infrared. The 
resulting datasets are first interpreted by comparing the observations with images and spectra of the close-by young cluster R136 
in the Large Magellanic Cloud and then by using model predictions for both the nebular emission lines and the spectral energy 
distribution of the sources. 

Results. We produce maps of the region containing the three sources in the R, J, Ks, L' bands and at 12.8 yum and perform their accurate 
relative positioning. We also provide spectra in the ranges 1. 8-2.4 yum, 3.3-4.0 //m, 8.1-9.3 fim and 10.4-13.2 fim. The spectral energy 
distribution of the three sources rises with wavelength. Emission lines from ionised hydrogen and molecular hydrogen are detected, 
as well as PAH emission. Conspicuous [Nell] 12.8 yum line is also present, while neither the [Arlll] 8.9 yum nor the [SIV] 10.4 yum lines 
are detected. This provides a stringent constraint on the age of the sources: we argue that they are relatively evolved young clusters 
(6-8 Myr). Owing to their ionising photon emission rates and ages, they must be extremely massive clusters (of the order of 10' Mq). 
Their mid-infrared spectral energy distribution suggests the presence of two components: (1) an optically thin component, with a 
continuum comparable to that of R136, and (2) an optically thick component which might be related to subsequent or on-going 
episodes of star formation. We anticipate that these sources are good candidates for evolving according to a bi-modal hydrodynamical 
regime, in which matter is trapped at the centre of a compact and massive cluster and generates further star formation. 

Key words. ISM: dust, extinction, ISM: HII regions. Galaxies: star clusters. Galaxies: individual: NGC1365, Infrared: galaxies 

1. Introduction little extinction. We prefer to use the generic term YMC in all 

cases, and specify, whenever needed, whether the YMC is still 



environments, suc h as dwarf ga laxies (O'Connell et al.'lg^i in Henize 2-1 (.Kobulnickv & Johnson 1 9991). in NGC5253 
Hunter etal. 1994: IO'Conndret al. 1995; Leitherer et al. 199'! (iGoriian et al.l l2()01h. in SBS 0335-052 (|Plante & Sauvagd 
Goriiani 1996i. interacting galaxies (Whitmore et al. 1993: 2002), in IIZw40 (Beck et al. 2002), in NGCl 365 and NGC 1808 



Goriiani 1996), interacting galaxies (Whitmore et al. 1993; 

Conti & Vacca.l994:,Shava et al... 1994; Whitmore & Schweizer ( Gallia ri o et al. 2(305; Galliano & Alloin 2008), in NGC7582 



Conti & Vacca.ib'^4- :^hava et al..,i ^^4;. Whitmore & :schweizer yjoiimnuc.^ ^uuj, v 
19951; lMeureijjT995l) and circumnuc l ear star-forming rings (|Wold & Galliano||2006|). 



('Benedict et al."l993t iBarth et al.lll995t iBower & Wilsonll 19951; As proposed in iJohnsonI (|2004 . parallel evolutionary se- 
[Maoz et al. 1996). Their possible parental link with classical and quences can be imagined for the formation of massive stars and 
well-studied globular clusters was proposed and, in subsequent for the formation of massive clusters. Both types of objects start 
studies, these young massive clusters (abbreviated YMCs) be- as HII regions deeply embedded in a dust cocoon. In the case 
came increasingly referred to as objects likely to evolve into of a massive star, this stage is identified as ultra compact HII re- 
globular clusters after a few Gyrs. Their masses are greater than gion (UCHII region) and in the case of a YMC as ultra dense HII 
10^ Mq, with radii smaller than 5pc and ages below 100 Myr region (UDHII region) (Kobulnickv & Johnson 1999). They are 
In the literature, the term super star cluster (SSC) is often asso- both inconspicuous in the visible and near-infrared (NIR), while 
ciated with YMCs which are bright in the visible, hence suffer bright in the mid-infrared (MIR) and far-infrared (FIR). They 
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are also intense sources of thermal radio continuum as well as of 
line emission from ionised gas. The embedding material even- 
tually dissipates off and they become detectable as UV-bright 
sources. The extent to which this parallel is sustainable remains 
an open question: similarities and differences along the two se- 
quences are worth investigating and might bring clues about the 
conditions for star formation in galaxies, and in particular for the 
formation of massive star clusters and globular clusters. 

We present in this paper the first step of a thorough analy- 
sis and modelling of embedded YMCs in nearby spiral galax- 
ies. The three YMCs we are interested in here are located 
in the starburst circumnuclear region of the Seyfert2 galaxy 
NGC1365 (dis t ance 18.6 Mpc, hence 1 " corresponds to 90pc). 
iGalliano et al.l (l2005l) discovered these sources in the MIR and 
demonstrated that they c oincide with bright the rmal centimetre 
radio sources detected bv lSandqvist et al. l (ll995l). They have also 
been r ecently detected in the CO molecule by ISakamoto et al.l 
(I2007h . These authors infer a mass of molecular material of the 
order of IO^Mq in the central 2kpc diameter region, which in- 
cludes thethreeYMCs under discussion. 

In iGalUano et al. (2005), their ages were grossly estimated 
from their radio spectral indexes and found to be of a few ( 3) 
Myrs. Accordingly, their masses were estimated to be of the or- 
der of IO'-Mq. 

These clusters are located at the inner Linblad resonance 
(ILR) in their host galaxy, and also within around 1 kpc from 
its active galactic nucleus (AGN). Their environment has a 
slightly above solar metallicity, as expected from the observed 
metallicity gradient in the disc of NGC1365 discussed by e.g. 
lDors&Copet tr(2005). 

The questions to be addressed in the current paper remain 
basic ones. What are the spectral characteristics of such YMCs? 
Are NIR and/or MIR data suflicient to get a trustable insight on 
the properties of embedded YMCs? Do we miss a substantial 
part of the phenomenon by observing only in the NIR? Which 
physical parameters can be safely derived for these objects, ei- 
ther in a direct way or through a comparison with models? 

The paper is structured as follows: the first part is devoted 
to the presentation and discussion of the data, while subsequent 
sections deal with the physical interpretation of the sources. In 
Sec. 12. II we describe the acquisition and reduction of the dataset, 
which consists of NIR/MIR images and spectra, all collected at 
the ESOA'LT using the instruments ISAAC and VISIR. We first 
discuss the images and perform the relative registration of the 
maps at different wavelengths (Sec. 12. 2t . We then describe and 
discuss the spectra (Sec. |2.3l l, which display intense nebular lines 
and NIR/MIR rising continua. In Sec. 12.41 we discuss the uncer- 
tainties on the quantities measured from the images and from the 
spectra. 

In Sec. [3] we perform a comparison of the data for the three 
YMCs in NGC1365 with comparable data obtained for R136 
in the LMC, one of the nearest known YMC. To do so, we 
use WFI and IRAC images, an ISO SWS spectrum and a wide 
field ISO CVF spectrum of R136. In SecH] we derive basic pa- 
rameters for the YMCs based upon their emission lines, via a 
comparison with predictions from a library of photo-ionisation 
models generated with the code CLOUDY. In Sec. |5] we analyse 
the NIR/MIR infrared emission of the sources, both their spec- 
tral energy distribution (SED) and their line emission, perform- 
ing a more complex modelling with the dusty stellar population 
evolution code GRASIL. Then, in Sec|6] we attempt to position 
our results in the light of the theoretical evolution of very mas- 
sive clusters, considering the bi-modal hydrodynamic sol ution 
for re-inserted matter, as proposed by ISilich et aTl (l2007h and 



iTenorio-Tagle et all ( l2007b . Finally, the conclusions and per- 
spectives of our work are highlighted in Sec. [T] 

2. The NIR/MIR dataset 

Using the ESO infrared facilities, we obtained new images 
of the central region of NGC1365 , as well as spectra of the 
three MIR/radio sources unveiled bv lGalliano et al.l ( l2005b . With 
ISAAC, the NIR spectro-imager at VLT/UTI, we have collected 
J (l.2idmE, Ks (2.2;um), L' (3.8;um) and M (4.5 yumfl images 
at an angular resolution of the order of 0.6 ", and low resolution 
long slit spectra of the three MIR/radio sources, in the K and L 
bands^. Notice tha t throughout this pap er, we retain the source 
nomenclature as in lGalhano et al.l ( l2005l) : the embedded sources 
are referred to as M4, M5 and M6. For one of the sources (M6), 
we also obtained a spectrum in the N band (around 10/vm]J with 
TIMMI2, the MIR spectro-imager of the 3.6m telescope at La 
Silla: despite its rather low S/N ratio in the continuum, a promi- 
nent [Nell] line could be detected. This encouraged us to per- 
form additional observations with VISIR at VLT/UT3, provid- 
ing an image in the narrow [Nell] filter at 12.8 //rrQ, together 
with low resolution spectra in the 8 fim, 1 1 fj.m and 12 fj.m bands 

B 

In addition, we use in our analysis a WFPC2 R band image 
retrieved from th e HST archive, the ATCA cen t imetr e maps by 
iForbes & NortsI ([1998) and by iMorg anti et a P ([1999). as well 
as the centimetre measurements of [S^andqvist e t al.l ( 19951) . We 
also consider the TIMMI2 10.4yL/m, 11.9yL/mand 12.9 //m images 
previously obtained bv .Galliano et al., (,2005i) . 



2.1. Data collection and reduction 

All data were acquired and reduced using standard techniques. 
Let us briefly recall the main steps below: 

ISAAC images: The images (J, Ks, L' and M bands) 
were obtained using the Aladin detector, with pixel scales of 
0.148" per pixel in the J and Ks bands and 0.071 " per pixel 
in the L' and M bands. The on-source integration times for the J, 
Ks, L and M bands were respectively of 90 sec, 300 sec, 300 sec, 
450 sec. A standard nodding technique was applied for the J 
and Ks observations, while chopping and nodding were used 
for the L and M observations. All images were reduced with 
the ECLIPSE package. The photometric calibration of the im- 
ages relied on observations of a standard star The precision on 
photometric measurements is of the order of 10%. The achieved 
angular resolutions are 0.56" in J, 0.44" in Ks, 0.39" in L' and 
0.38 " in M. The M band image has a low S/N ratio though and 
is not used in the following. 

ISAAC spectra: The low resolution long slit mode was used 
to collect Ks and L band spectra of the sources M4, M5 and 
M6. Spectral resolutions were respectively R=450 and 360 for 
the Ks and L bands. Two positions of the 1 " width slit al- 
lowed to obtain the three spectra: one slit at PA=2.8° (PA pos- 
itive from North to East) passed through the AGN and M4, 
while another sHt at PA=145.r passed through M5 and M6. 
T he shts were pr ecisely positioned following the measurements 
in IGalliano et al.l ([2005.) . and by performing blind offsets refer- 
enced on the AGN. Standard data reduction procedures were ap- 



' P074.B-0166; October 29-30, 2004 
- P072.B-0397; December 01-02, 2003 
3 P074.B-0166; October 23-24, 2004 
P074.A-9016; December 01, 2004 
^ P076.B-0374(A); November 18 & 20, 2005 
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plied, using ECLIPSE and IRAF. In the Ks band, self-chopping 
effects prevent assessing precisely the continuum level, although 
they do not affect the measurement of emission line fluxes. The 
spectra were extracted through slit windows of 1.4" along the 
slits. Night-sky lines were used for the wavelength calibration. 

TIMMI2 spectrum: With TIMMI2, an N band spectrum 
was obtained for M6, using the lOjum low resolution grism 
(7.5 //m to 13.9;um) with spectral resolution R=160. The pixel 
scale is 0.45". The slit, 1.2" wide and at PA= 0°, was blindly 
positioned on the target using the AGN as a reference. The stan- 
dard nodding and chopping technique allowed an efficient re- 
moval of the background emission. The on-source exposure time 
was about 4000 sec. Subsequent data processing consists in the 
addition of chopping and nodding pairs, followed by a shift-and- 
add procedure to sum the two negative and one positive spectra 
generated by the chopping-nodding technique. 

VISIR iitmgc: VISIR allowed to collect an image of 
NGC1365 through the narrow [Nell] 12.8 yum filter (FWHM of 
0.2 /im). The standard chopping and nodding technique enabled 
to remove the MIR background. The data reduction consists in 
shifting and stacking the individual frames (each corresponding 
to a chopping position). The achieved angular resolution is 0.4 ". 
To optimise the detection of extended features, the image was fil- 
tered using th e mr .filter routine of th e MR/1 software package 
developed by Murtagh & Starckl (Il999l) . 

VISIR spectra: We chose for the VISIR spectra the same 
two slit positions as for the ISAAC spectra. We observed in 
three settings with centred wavelengths at 8.5/vm, 9.8 jum and 
12.2 jum , in a low resolution mode (/l/A/l of the order of 200). 
Exposure times were respectively 655 sec, 510 sec and 516 sec 
for the three settings. Doing so, we achieve the following spec- 
tral coverage: from 8.1jum to 9.3 //m, and from 10.4 /^m to 
13.0jum, with a good overlap around lljum. Observing pro- 
cedures and data reduction techniques are similar to those de- 
scribed for the TIMMI2 spectra. The absolute calibration of the 
12.2jum centred spectrum was assessed through the calibrated 
VISIR 12.8 /im narrow band image. The 9.8jt/m centred spec- 
trum was then scaled, through its overlap with the 12.2 /im cen- 
tred spectrum. For the spectrum centred at 8.5 //m , we retain 
the original flux calibration. Therefore, there is some uncer- 
tainty between the relative fluxes of the 8.1-9.3 yum and the 10.4- 
13 fim segments of the spectrum. The precision on the absolute 
calibration of the 10.4 to 13 fim segment is of the order of 20%, 
comparable to that of the VISIR image and better than that of 
the 8.1-9.3 fim segment. The calibration uncertainty for the lat- 
ter spectral segment is mostly due to possible varying slit losses, 
and might reach 40%. 

2.2. Analysis of the images and registration 

Fig.[T]displays a three colour image (0.6 //m, 2.2 //m and 3.8 yum) 
of the central region of NGC1365. The northern dust lane ap- 
pears in green (2.2 fj.m), while the three sources of interest show 
up in red. Fig. |2] provides a summary overview of the imag- 
ing data. The large image presents the central 30"x30" of 
NGC1365 through the HST F606W filter (referred to as an R im- 
age since it is close to the standard R band filter). This image is 
the pipeline-reduced HST/ WFPC2 archive image. The location 
of the bright type 2 AGN ( LindblaJ 19991) is at relative coordi- 
nates (0 "; "). The square to the North of the AGN outlines the 
region where the three bright MIR/radio sources, M4, M5 and 
M6, are found. The series of small images offers closer views of 
this region at wavelengths: 0.6 jum (WFPC2), 2.2 fim (ISAAC), 
3.8 yum (ISAAC) and 12.8 fim (VISIR), as well as a reprint of the 




Fig. 1. Three colour image of the central region of NGC1365. 
The WFPC2 R band image is blue, the ISAAC Ks band image is 
green, and the ISAAC L' band image is red. North is up and East 
is to the left. The dimensions of the image, 25 "x30 "correspond 
to 2.25 kpcx2.7kpc for a distance of 18.6Mpc 




Fig. 4. Comparison between the VISIR narrow b and 12.8 //mfil- 
ter image (left) and the 3.5cm map (right) by iMorganti et al.l 
( 1999) for a 10 "x 14 "region including the AGN (bottom source 
on the left image), the three embedded YMCs and the MIR/cm 
ring. North is up. East is to the left. The ticks along the axis are 
drawn with 0.5 " intervals. 



3 cm ATCA image from lForbes & Norrisl(ll998l) . On these closer 
views, three pairs of concentric circles pinpoint the positions of 
M4, M5 and M6, and the aperture sizes used for the flux density 
measurements presented in Sect. 12.41 

The morphology of the nuclear region of NGC1365 varies 
with wavelength. Hence, the relative registration of images at 
different wavelengths has to be performed with care. In the fol- 
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Fig.2. The embedded clusters in NGC1365. The large image shows the inner 30"x 30" observed with HST/WFPC2 in the 
F606W filter (0.6 /vm). The highlighted square shows the region where the three embedded young massive star clusters M4, 
M5 and M6 are located. The small images detail this regio n at five waveleng ths: 0.6 /urn (HST/WFPC2), 2.2 fim (VLT/ISAAC), 
3.8 jum (VLT/ISAAC), 12.8 jum (VLT/VISIR) and 3 cm ( Forbes & Norrislll99i ATCA, reproduced from their work). For clarity, 
the J (1.2 fim) image mentioned in the text is not displayed here. Its angular resolution is rather low and adding it would not bring 
any pertinent information. The locations of M4, M5 and M6 are given by the centres of the concentric circles on the small images. 
They correspond to the locations of the 12.8 ;um sources. Sources labelled CI, C2 and C3 as well as the active galactic nucleus 
(labelled AGN) are used for the registration (see Sect. 12.2b . 



lowing, we do not provide absolute source positions but rather 
derive relative registrations. For image orientations and pixel 
sizes, we use the values attached to the image headers. On the R, 
J and Ks band images, enough sources are detected in the field 
of view to allow good relative positioning. For this, we use the 
sources labelled CI, C2 and C3 in Fig.|2]which lie outside zones 



of highest extinction. We do not use the AGN itself to register 
these images, as a shift between its visible and infrared peak- 
positions may occur. On the other hand, on the L' and N band 
images, only the AGN and the sources M4, M5 and M6 have a 
S/N ratio sufficient to perform positional measurements. In this 
case, the AGN must be used for registration and we make it co- 
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L' 12.8/xm 



CI 



C2 




1.4 1.2 1.0 o.e 0. 



C3 




7.6 7.4 7.3 7.0 6.E 



AGN 




□□□□ 



0,4 0,2 0,0 -0,2 -0,-^ 



M4 



M5 




M6 

Fig. 3. Details of the different images illustrating, for CI, C2 and C3 the quality of the registration, and for M4, M5 and M6, the 
morphological changes that these sources suffer with wavelength. Each raw shows the images of a given source at the various 
observed wavelengths. The wavelengths are from left to right: R, J, Ks, L' and 12,8 fim. The sources are from top to bottom: CI, 
C2, C3 (see Fig,|2]i, AGN, M4, M5, M6, The first four raws illustrate the quality of the registration. The last three raws highlight the 
morphological changes of the embedded YMCs with wavelength. The X and Y scales for each image represent the same coordinates 
as on Fig. 12] where (0"; 0") coincides with the AGN location. 



6 



Galliano et al.: Extremely massive star clusters in NGC1365 





Fig. 5. Frames (a, b, c) from left to right: R, Ks and 12.8 fim images of the embedded clusters region. North is up and East is to the 
left. The slit drawings show the positions of slit 1 (through the nucleus and M4) and slit 2 (through M5 and M6). Along the slits, 
the thin and short ticks mark a scale in arcsec, while the thick and long ticks mark the position where the spectra are recorded. The 
references are the nucleus for slit 1, and M5 for slit 2. Frames (d & e): ISAAC spectra through slit 2 after reduction and before 
extraction. The wavelength increases from left to right. In the Ks spectrum, the two brightest lines are Pa and Bry. This image 
shows that, due to self-chopping, it is difficult to assess the actual continuum level. In L', the bright line to the right is Bra. Since, 
no extra emission is detected along the slit -apart from the two clusters- self chopping is inexistent. Frame (f) VISIR spectra, from 
left to right separated by grey bands: 8.5 //m setting (8. 1-9.3 //m), 9.8 jum setting (10. 4-12.4 yum), 12.2, yum setting (11.3-13.0jum). 
The bright emission line on the right is [Nell]. 



incide with its location on the Ks band image (Fig. [3]l. On the 
grounds of modelling, this is an acceptable assumption as only 
a minor offset, if any, is to be expected for the AGN position 
between the Ks and N bands dGranato et alj|1997l) . 

Under such a registration, Fig.[3]ascertains that the positions 
of the reference sources CI, C2 and C3 agree within better than 
+0.1 " on the R, J and Ks band images. At longer wavelengths, 
as discussed above, one must rely on the location of the AGN 
and, assuming that offsets of the nuclear peak projected posi- 
tions in the Ks, L' and 12.8/vm bands are small compared to 
the image resolutions, the global precision on the registration 
of all maps is better than ±0. 1 ". Estimate of the precision in the 
relative positioning does matter since, as shown on the three bot- 
tom rows of Fig. [3] the positions of the peaks corresponding to 
M4, M5 and M6 appear to move slightly with wavelength. In the 



following, we call "location" of M4, M5 or M6 the projected po- 
sition of their corresponding 12.8 fj.m peak, which are provided 
in table[T] They correspond to projected distances to the AGN of 
640 pc, 920 pc and 760 pc for M4, M5 and M6 respectively. 

In the R and Ks bands, no conspicuous emission is detected 
at the exact location of M4. Notice that M4 is located close 
to, but not coincident with, the apex of a cone-shaped struc- 
ture appearing on the R band image, roughly extended along the 
North-East direction. This structure is also bright in the J band. 
Emission in the Ks band is detected between the location of M4 
and the peak in the J band mentioned above. It may be a mixture 
of emission from M4 and from the "cone". The source M4 is 
clearly detected in the L' band as well as at 12.8 /urn. 

The source M5 is detected in the R, J, Ks and L' bands and 
at 12.8 jum. 
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M4 M5 M6 

value value value 

low high low high low high 
position with respect to AGN [arcsec] 



Aa [arcsec] 


0.31 


0.51 


2.70 


2.90 


4.70 


4.90 


A(5 [arcsec] 


7.02 


7.22 


9.76 


9.96 


6.80 


7.00 






flux densities [mJy] 






R 


<0.10 


0.13 


0.19 


<0.05 


J 


<1.21 


0.79 


1.88 


<1.01 


Ks 


<1.18 


0.62 


1.47 


0.31 


1.29 


L' 


0.28 


3.80 


0.69 


3.87 


1.81 


7.16 


10.4 yum 


0.57 


3.94 


41.33 


62.90 


37.18 


57.01 


11.9yum 


20.10 


47.26 


78.76 


130.16 


109.15 


169.90 


12.9 yum 


91.97 


169.03 


208.83 


341.28 


379.61 


575.43 


12.8yum 


77.24 


128.32 


183.10 


287.64 


264.64 


421.92 


6cm 


2.38 


2.98 


1.12 


1.72 


2.45 


3.05 


20cm 


3.04 


3.64 


2.44 


3.04 


2.24 


2.84 






line fluxes [10 erg s 


' cm^] 






Pff 


23.18 


34.76 


64.64 


96.96 


60.43 


90.65 


Bi6 


1.23 


1.84 


3.40 


5.10 


3.48 


5.22 


H. 1-0S(3) 


0.71 


1.07 


0.96 


1.44 


1.29 


1.94 


H, 1-0S(2) 


0.36 


0.54 


0.37 


0.56 


0.44 


0.66 


Hel 


0.86 


1.29 


1.85 


2.77 


2.44 


3.67 


H2 l-OS(l) 


1.06 


1.59 


1.02 


1.53 


1.13 


1.70 


Bry 


2.71 


4.06 


5.87 


8.81 


6.84 


10.26 


H2 l-OS(O) 


0.38 


0.57 


0.43 


0.64 


0.52 


0.78 


H, 2-lS(l) 


0.27 


0.40 


0.31 


0.47 


0.31 


0.46 


H, l-OQ(l) 


0.40 


0.60 


0.74 


1.10 


1.27 


1.91 


H, 1-0Q(2) 


0.33 


0.50 


0.44 


0.66 


0.52 


0.78 


H, 1-0Q(3) 


0.37 


0.55 


0.52 


0.78 


0.58 


0.87 


m 


38.70 


58.06 


35.76 


53.64 


59.23 


88.85 


Bm 


16.66 


24.98 


19.82 


29.74 


31.50 


47.24 


[Arlll] 




<12 


<11 


<12 


[SIV] 




<6 


<5 






<5 


[Nell] 


110 


170 


440 


660 


570 


850 






"spectral indexes Fv oc v " 








-0.49 


-0.29 


-0.87 


-0.67 


-0.54 


-0.34 


-0.29 


-0.09 


-0.66 


-0.46 


-0.03 


0.17 



" spectral indexes from lSandqvist et^ ( Il995h . An uncertainty of ±0. 1 is added. 



The source M6 is detected in the Ks, L' bands, as well as at 
12.8 fim. Possible counterparts to M6 are detected, slightly offset 
to the South, in the R band (offset by 0.2 ") and J band (offset by 
0.1"). 

Regard ing the 3 cm image, we o nly h ave access to the printe d 
figures bv iForbes & NorrisI (Il998l) and iMorganti et al] (1 19991) . 
therefore no precise measurement can be performed. Fig. |4] 
presents the filtered 12.8 fim image with cuts highlighting the ex- 
tended emission and compares this image to the 3.5 cm map by 
iMorganti et al.l (Il999b . The similarity of the two maps is striking, 
not only for the main peaked sources, but also for the extended 
emission and in particular for the western feature elongated in 
the North-East direction. The sources M4, M5 and M6 appear 
to be distributed in a well defined star-forming ring, delineated 
through its MIR/radio emission. The AGN itself is neither radio 
bright, nor clearly isolated on the radio image, and hence cannot 
be used for regis tration purp oses. The two independ ent maps by 
iForbes & NorrisI ([T998i) and IMorganti et~an d 19991) are in good 
agreement. Both indicate that the three radio sources cannot be 
made precisely coincident with the three MIR sources simulta- 
neously. In Fig.|2l M6 is chosen as the reference: M4 and M5 are 
slightly shifted to the North-East with respect to the radio map. 
In Fig.m we display an alternative solution in which M4 and M5 



have precise radio counterparts: then, M6 appears to be slightly 
shifted to the South- West with respect to the radio emission. As 
each of the two independent radio maps shows this offset -in its 
comparison with the MIR image-, one is led to conclude that the 
offset is real, but remains to be understood. 

The sources M4 and M6 appear to be located on the dust 
lane, in projection. The fact that, in the R band, M4 is not de- 
tected and M6 is very weak suggests that they lie inside the dust 
lane. The source M5 lies on the edge of the dust lane, and con- 
sistently appears less extincted, showing a bright counterpart in 
the R band. Halos of emission are detected around M4, M5 and 
M6, both in the L' band and at 12.8 yum. 

In conclusion, the relative positioning of the visible and IR 
images has been achieved at a precision of 0.1 "and is adopted 
in the following analysis. 

Let us make some cautionary remarks in the case of M4: the 
comparisons between its maps indicate that aperture flux mea- 
surements should be unreliable in the R, J and Ks bands, where 
no clear counterparts are found and as it lies in the vicinity of an 
R band emitting cone and of a J to Ks band emitting region. We 
believe that, at these wavelengths (R, J and Ks bands), aperture 
flux measurements for M4 are upper limits. 
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Fig. 6. From top to bottom, Ks, L' and N band spectra of M4. Fig. 7. From top to bottom, Ks, L' and N band spectra of M5. 



The source M5 is well defined at all wavelengths. A neigh- 
bouring weak source to the North (seen on the R band image) 
may contribute some flux to its aperture flux measurements in 
the J and Ks bands. Yet, M5 probably remains the dominant 
source in the Ks band. 

The source M6 is clearly detected at all wavelengths except 
in the R band. At this wavelength, a weak source is detected 
-0.2" to the South of the location of M6 (See Fig.|2]and Fig.EJ. 
We cannot firmly exclude that it is an unrelated source, therefore 
its flux has been included in computing the R band flux upper 
limit of M6. 

These remarks of course apply to continuum measurements 
from the spectra that we are about to discuss, since they corre- 
spond to a 1 "xl.4" aperture. This leads to the following con- 
clusion, to be kept in mind when interpreting the NIR spectra: 
in the case of M4, the Ks band continuum is an upper limit, in 
the case of M5 it may be slightly overestimated, and in the case 



of M6 it is a proper estimate. At longer wavelengths, such issues 
are irrelevant. 



2.3. Analysis of the spectra 

In the upper part of Fig.|5] we provide drawings of the slits used 
for the spectroscopy in the Ks and L' bands, superimposed on 
images in the R and Ks bands and on the 12.8 jum image. Along 
the slits, thin ticks are drawn every arcsec, while longer and 
thicker ticks highlight the positions at which emission line spec- 
tra have been recorded. References used for the positioning of 
the spectra along the slits are the AGN for slit 1 and M6 for slit 
2. M6 is a good reference since its position is well defined and 
does not shift significantly with wavelength above 1 fim. 

For M4, the emission line spectrum position (given by the 
thick tick mark on the slit drawing) coincides with the emission 
peak on the 12.8 //m image; it is located 0.1-0.2" to the South of 
an emission knot in the Ks band (but has itself no clear counter- 
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Fig. 8. From top to bottom, Ks, L' and N band spectra of M6. 



part in Ks); it is also located 0.2 " to the East of the apex of the 
cone-shaped structure seen in the R band image. 

The emission line spectrum at the position of M5 corre- 
sponds to well defined emission knots on the three images (R, 
Ks and 12.8 /zm). 

The emission line spectrum at the position of M6 corre- 
sponds to an emission knot both on the 12.8 fim and Ks images, 
but is 0. 1-0.2 " to the North of an emission knot in the R image. 

The 2D raw spectra obtained through slit 2 in the Ks and L' 
bands with ISAAC, and in the N band with VISIR, are displayed 
at the bottom of Fig. |5] The brightest emission lines identified 
on the figure are Pa and Bry in the Ks band. Bra in the L' band 
and [Nell] in the N band. 

The extracted spectra are displayed in Figs. |6] [7] and [8] for 
M4, M5 and M6 respectively. The following set of emission 
lines is detected: (a) nebular lines Pa, Bt6, Hel, Bry, Pf6, Pfy , 
Brff, [Nell] 12.8 yum, (b) molecular lines from H2 and (c) PAH 
emission. We detect PAH signatures at 3.3 jum and 11. 3 fim. The 
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Fig. 9. Position- Velocity {Pa line) diagram for M4 along slit 1 
(PA=-2.8°). The X-axis represents the velocity, the X-0 value 
being set at the emission peak. The Y-axis represents the spatial 
direction, top corresponding to the North. The Y-0 coordinate 
corresponds to the emission peak, also marked by a thick tick 
along the drawing of slit 1 in Fig.|5] At any position, the instru- 
mental broadening is dominant (660 kms"'). 



bump-like feature on the 8-9 jum spectrum of M6 might be a sig- 
nature of the 8.6 //m PAH, although it looks a bit broad. Even 
though the [Nell] 12.8 ;um line is bright, we do not detect any 
[ArIII]8.9;um or [SIV]10.5yL(m line emission, which is a puz- 
zling and interesting result. The 12.7 yum PAH feature, expected 
to show up as a broad feature, is also absent. For the three 
sources, the continuum spectral distribution is flat or slightly ris- 
ing with wavelength in the Ks band and neatly rising in the L' 
and N bands. The presence of the 9.7 fim silicate absorption band 
is inferred in M4, M5 and M6, from the simultaneous fall and 
rise of the 8.1-9.3//m and 10.4-13. 0/im continua respectively 
matching the blue and red wings of the silicate absorption fea- 
ture. For the continuum in the Ks band, one must bear in mind 
the remarks given at the end of the previous section. Between 
3 fim and 3.3 yum , the transmission of the atmosphere is poor: 
this part of the spectrum is noisier and we think that the uncer- 
tainty on the atmospheric correction induces a supplementary 
error on the continuum value. We are tempted to believe that 
only the segment of the spectrum red-wards of 3.3 /im is reli- 
able. In the N band, the spectral calibration for the red side of 
the spectrum has been performed through scaling to the VISIR 
narrow band 12.8 fim image and taking advantage of the overlap- 
ping region between 1 1 fim and 12 fim. This procedure resulted 
in applying scaling factors to the spectra. Such a correction can- 
not be applied to the blue side of the N band spectrum, since no 
overlapping spectrum is available. Therefore, the relative scaling 
between the 8.1-9.3 fim segment and the 10.4-13 fim segment of 
the spectrum may be uncertain by a factor up to 2. 

On the spectra, an interesting feature is the extension of the 
line emission along the slit direction (see the Ks band spectrum 
in Fig.|5]l for the brightest lines. For M4, and in spite of a rather 
low spectral resolution, we can even detect a velocity gradient in 
the Pa line. The position-velocity diagram for the M4 Pa emis- 
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Fig. 10. NIR/MIR SEDs of the embedded clusters M4, M5 and 
M6 in NGC1365 (red points with error bars or upper limits), 
over-plotted on the Ks, L' and N band spectra (black continu- 
ous line). The NIR me asurements are from t his article, while the 
MIR points are from iGalhano et all (l2005h . For M4, since no 
direct Ks band counterpart is observed on the image, the photo- 
metric point is given as an upper limit. This apparent inconsis- 
tency denotes the fact that we do not formally associate the Ks 
band continuum to the cluster emission (see text). 



sion line is displayed in Fig.|9l Between positions -1 " and +1 ", 
a velocity difference of the order of 100-200 km s"' is measured. 
This kinematical feature will be discussed elsewhere. In the case 
of M5 and M6, we also detect extended nebular emission (see 
Fig.|5]l, but no velocity gradient. 

Spectral analysis in the radio d omain would also be of great 
interest. The radio measurements bv lSandqvist et al.l (1 19951) indi- 
cate that the spectral indexes in the three sources flatten with in- 
creasing wavelength, suggesting optically thick radio emission. 
Unfortunately, such measurements are delicate from an image 
wi th complex brightness distribution and have not been repeated 
bv lMorganti et al.1 (Il999l) . hence cannot be double checked. 

In summary, the collected spectra exhibit bright nebular 
and molecular emission lines at the location of the MIR/radio 
sources M4, M5 and M6: this reinforces the idea that these 
sources are embedded YMCs. The line emission is spatially ex- 



tended; in the case of M4, a velocity gradient over a few arcsec 
is detected, suggesting the presence of an outflow. 

2.4. Measurements and related uncertainties 

Table [T] displays the full set of measurements performed for 
the three sources M4, M5 and M6. The flux density measure- 
ments were made through apertures of 0.6 " radius, represented 
in Fig.|2]by the small inner circles, and centred on the source po- 
sitions as determined on the VISIR 12.8 fim image. Fluxes were 
measured on the following images: the HST/WFPC2 R band im- 
age, the ISAAC J, Ks, L' band images, the TIMMI2 10.4jum, 



ri.9yum and 12.9 /zm im ages ( Galliano et alj 20051) . and the 
VISIR 12.8 fim image. In iGaUiano et alj (|2005I)T we suspected 



the 12.9 yum flux densities to be overestimated by a factor 2. The 
TIMMI2 N band spectrum allows a direct measurement of the 
1 1.9 yum/12.9 jum flux ratio (independent of the flux calibration), 
and shows that our suspicion was found ed indeed. Therefore , 
die 12.9 fim fluxes for NGC1365 quoted in lGaUiano etal] (l2005h 
must be decreased by a factor 2. 

Around each aperture, a "background" level was estimated 
by computing the median pixel value in an annulus with radii 
0.6" and 1.2" around the aperture centre. The large and outer 
circles on Fig. |2]represent the outer borders of such annuli. We 
cannot be certain that this "background" also affects the source 
within the small aperture, as its origin and its location with re- 
spect to the source are unknown. Therefore, in a conservative 
approach we provide two flux density measurements for the 
sources M4, M5 and M6: one after "background" subtraction 
and one without "background" subtraction. The two figures for 
each flux represent the measurement uncertainty resulting from 
an ambiguity in the interpretation. To this uncertainty we add 
quadratically the uncertainty due to the photometric calibration. 
The calibration uncertainties used for the different bands are as 
follows, R: 10%, J: 10%, Ks: 20%, L': 20%, 10.4 //m: 20%, 
11.9yum: 20%, 12.9jum: 20% and 12.8 A^m: 20%. Whenever at 
a given wavelength no counterpart is clearly identified at the po- 
sition of the MIR source, then we consider the high value of the 
error bar as an upper limit. In table [1] the low and high values 
given for each measurement correspond to the lower and upper 
limits for the given measurement due to the uncertainties defined 
above. 

Table [1] also provides the flux measurements for the emis- 
sion lines identified in the spectra of M4, M5 and M6. In the 
case of M6, we have additional measurements for lines in the 
N band. Fina lly, we recall in t his tab le the radio measurements 
published by ISandqvist et alj (Il995h . Fig. [10] displays, for the 
three embedded YMCs, their SED together with their combined 
spectrum, covering in total a range from 0.5 jum to 13/im. Let 
us briefly discuss the comparison between SEDs and spectra. 
The SED points show flux density values that we believe cor- 
respond directly to the sources of interest. However, in some 
cases they may correspond to upper limits because of source 
confusion. On the contrary, the spectra contain effective signal 
from the source and its surrounding. This distinction is partic- 
ularly relevant in the case of M4. That is why on Fig. [TOl for 
M4 the Ks band spectrum is shown as an upper limit. For the 
sources M5 and M6, the Ks band continuum corresponds to the 
upper part of the error bar on the flux measurement, because 
this value corresponds to the measurement performed without 
"background" subtraction. The L' band data points and the N 
band data points of Galliano et al. (2005) are in excellent agree- 
ment with the spectra. The TIMMI2 10.4 jum flux density for M4 
is not reported on this figure, since there was no detection on the 
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TIMMI2 image and the, then quoted, upper limit of 5mJy may 
have been underestimated. 



3. Comparison with a resolved young massive 
cluster: R136 in the LMC 
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Fig. 11. Extinction curve used in this paper This curve has 
been derived with G RASIL for the average Galactic extinction 
dPraine & Leelll984 . 



1000.000 



100.000 



c 

T3 



10.000 



1.000 



0.100 



0.010 



0.001 




6 8 10 

wavelength [/im] 

Fig. 13. Comparison of the ISAAC/VISIR spectrum of 
NGC1365/M6, the ISO CVF spectrum of R 136 and ISO SWS 
spectrum of R136. The spectra of R136 are projected to the dis- 
tance ofNGC1365. 



The fact that the MIR/radio sources M4, M5 and M6 in 
NGC1365 exhibit nebular lines strongly supports an interpreta- 
tion in terms of embedded YMCs still hosting a fair fraction of 
ionising stars, hence younger than lOMyr. 

Let us start the analysis by comparing the available data for 
the YMCs in NGC1365 with similar information (images and 
spectra) on R136, a young cluster in the LMC. Basic parameters 
for these clusters will be derived in subsequent sections. 

R136 is one of the nearest example of a young massive 
cluster It lies in the most luminous star forming region in the 
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Fig. 15. Predicted line fluxes for the model described in Sec.E 
From top to bottom: Bry, Bra, [Arlll], [SIV] and [Nell]. The X- 
axis displays the age in Myr and the Y-axis provides the line flux, 
at the distance of NGC1365, corresponding to a IQ^Mq cluster 



Local Group: 30 Doradus in the LMC, at a distance of 50kpc. 
R136 has an age in the range 1-3 Myr and a stellar mass of 
SxlO'^Mo (Boulanger & Rubio 2006). 

Even though R136 is not properly speaking an embed- 
ded cluster, it is an interesting object for comparing with the 
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Fig, 12. ESO 2.2m WFI/Spitzer IRAC view of R136 in the Large Magellanic Cloud. The field of view is lO'xlO'. North is up and 
East is to the left. This region (145pcxl45pc) would project onto a field of 1.6 "xl.6" at the distance of NGC1365. Blue codes the 
WFI Ha image, green the WFI V image and red the IRAC S.Oyum Spitzer image. The 3' scale corresponds to the aperture size of 
the CVF spectrum aperture and the arrow shows the location of the 14"x20" ISO SWS spectrum aperture. 



MIR/radio emitting sources of NGC 1 365 . A three colour map of 
R136 is displayed in Fig. [12] built as a composite of two archived 
images obtained with the Wide Field Imager (WFI) at the 2.2m 
telescope on La Silla (ESO) and one archived image collected 
with IRAC on board the Spitzer space telescope (NASA). Green 
codes the WFI V-band image (tracer of the stellar component), 
blue codes the WFI Ho- image (tracer of the ionised gas) and 
red codes the IRAC S.Ojt/m image (tracer of dust). The resolu- 
tion of the WFI images has been degraded to 3.2", which is the 
resolution of the IRAC 8.0/^m image. The composite field size 
is lO'xlO', which corresponds at the distance of the LMC to 
145xl45pc^. This field projects into a region of 1.6 "xl.6" at 
the distance of NGC 1365. This image shows that, even if R 136 
does not qualify as "embedded" cluster, it is a bright source of 
MIR radiation: would it lie behind or within an extinction lane, it 
would indeed be identified as an "embedded" cluster This sug- 
gests that the MIR emitting dust of a so-called "embedded" clus- 
ter may not be necessarily in direct relation with the "embedded- 
ness" of the cluster, as assumed in the typical cocoon-like view 
of an embedded cluster. As highlighted on this colour image, 
the three components are spatially separated: one can clearly 
identify the Ha emitting gas, spatially correlated with but dis- 
tinct from the MIR dust component, and the stellar cluster itself. 
The wide dust shell (bubble) pushed by the cluster is prominent, 
while the Hor emission shows up as a thin layer at its inner rim. 



Of the many spectra available for R136, let us consider: (i) 
the ISO SWS spectrum of a 14"x20" region in the North-East 
MIR emitting lobe (shown by an arrow on Fig. [T2T i and (ii) the 
ISO CVF spectrum of the 3'x3' region around the central clus- 
ter, covering the wavelength range from 2 to 16 /vm (taken from 
Fig. 3 of .Boulanger & Rubio 2006) . 

Notice that an aperture of 3'x3' in the LMC translates 
to 0.5"x0.5" at the distance of NGC1365. This allows a di- 
rect comparison between the CVF spectrum of R136 and the 
ISAAC/VISIR spectra of the embedded YMCs in NGC 1365. 
Fig.[T3]displays the SWS and CVF spectra of R136 projected at 
the distance of NGC1365, as well as the spectrum of N1365/M6. 
This comparison shows that M6 is, intrinsically, almost two or- 
ders of magnitude brighter than R136 and that, given the avail- 
able data, the two sources exhibit at first order a rather simi- 
lar spectrum. A remarkable diff'erence, though, is the fact that 
in the 8-13 //m wavelength range, three bright emission lines 
are detected on the SWS spectrum of R 136 ([Arlll] at 8.9 //m, 
[SIV] at 10.5 /im and [Nell] at 12.8 yum) while only the [Nell] 
12.8 fim line is observed in the case of M6. This diff'erence holds 
true in the case of M4 and M5. Owing to its low spectral reso- 
lution the R136 CVF spectrum does not show any of the narrow 
forbidden lines but only broad features such as the PAH bands. 

We provide in Fig. [141 ^ comparison of the WFPC2 R and 
the [Nell] 12.8 yum images of M5 with the WFI Ha and the IRAC 
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Fig. 14. Map comparison between NGC1365/M5 (left column) and R136 projected at the distance of NGC1365, The R136 WFI 
Ha image and IRAC 8 jjm image, after projection at the distance of NGC1365, have been degraded to the seeing and pixel size of 
the NGC1365/M5 WFPC2 R image and the VISIR 12,8 fj.m image respectively. 



8,0jum images of R136 projected at the distance of NGC1365, 
with same seeing and same pixel size as the M5 images. For this 
comparison we have chosen M5 since, among the three YMCs 
in NGC1365, it is the one with the brightest visible counter- 
part, hence the closest to the evolutionary stage of R136, On 
the left side of Fig, [141 the images of M5 are displayed, while 
on the right side, the corresponding degraded images of R136 
are shown. Notice that the bulk of the 8,0 yum emission in R 136 
looks elongated, while its counterpart in M5 is unresolved. The 
same occurs with the visible image: at the resolution of HST, 
the degraded R136 looks slightly more extended than M5, This 
comparison suggests that, in terms of intrinsic properties, the 
YMCs in NGC1365 are sHghdy smaller than R136 while they 
are about two orders of magnitude brighter: they definitely de- 
serve to be called "compact". Another piece of evidence comes 
from the HST image, where M5 remains unresolved, at a reso- 
lution of 0, 11 ", which corresponds to about lOpc, The fact that 
the sizes of the YMCs in NGC1365 are comparable to, or even 
less than, that of R136, which is 100 times fainte r, supports th e 
idea that star cluster sizes do not depend on mass (lLarsenll2004l) . 

This simple and direct data comparison between the 
MIR/radio sources in NGC1365 and R136 adds support to the 
interpretation of the NGC1365 sources in terms of compact 
YMCs. 

An additional striking difference between the two sources is 
the absence of the [Arlll] and [SIV] line emission in the YMCs 



in NGC1365. This important difference is discussed in Sec. |4] 
and explained through an age effect, the clusters in NGC1365 
being older than R136. A detailed comparison of the SEDs of 
R136 and of the YMCs in NGC1365 is performed and included 
in Sec.lH 



4. Derivation of the YIVIC basic parameters 

In this section, we attempt to derive parameters for the YMCs in 
NGC1365, such as extinction, age and mass, using their emis- 
sion line flux measurements. We make the gross assumption that 
a uniform foreground emission is affecting in a similar way all 
the emitting components in the cluster This is obviously an over- 
simplification, but it leads to first-order interesting conclusions 
while avoiding to use model-fitting with too many free parame- 
ters. We shall turn to a more realistic modelling in the next sec- 
tion. 

4.1. Extinction 

First, let us use the Bry and Brff line measurements to derive 
the extinction towards the YMCs. We know already that they 
are located in a region of large extinction, in the vicinity of, or 
within the prominent East- West dust lane to the North of the 
AGN (Fig.|2]l. In the following discussion, we use the extinction 




Fig. 16. For each cluster, and for the extinction derived in Sec. 14.11 this figure gives the ratios between observed and modelled 
[Arlll], [SIV] and [Nell] lines after correcting for the extinction correction and normalising the modelled Bry fluxes to the observed 
ones. The ratio is set to one when the modelled flux is inside the error bar, or below the upper limit, which is the case for [Arlll] and 
[SIV]. The total stellar mass (taking into account evolutionary eff'ects) is also plotted for each possible age. 



curve derived with GRASIL ("Sil va et al.lll998l) for Galactic dust. 
This curve is shown in the NIR/MIR range in Fig. II l1 

A direct estimate of the extinction towards the embedded 
cluster nebular gas can be retrieved from the Bry/Bro- line ra- 
tio, assuming optically thin conditions for the nebular gas itself. 
This ratio depends only mildly on the gas density. Moreover, 
considering that the temperature of HII regions tends to cluster 
around Tg = 10"* K (Osterbrock 1989), the predicted line ratio of 
Pa ( 1 . 87 fim) to Bt6 ( 1 .94 fim) is 0.0545 . The diff'erential redden- 
ing between Pa and BiS is small, the two lines being quite close 



in wavelength, so we can use the directly measured PafBrd ratio 
to check the conditions in the gas, and the validity of assuming 
the case-B line ratios. For the embedded clusters M4, M5 and 
M6, the measured ratios are respectively 0.053, 0.053 and 0.058, 
in good agreement with the theoretical value 0.0545 mentioned 
above. 

Consequently, the case-B assumption looks adequate and we 
can confidently derive the extinction towards the embedded clus- 
ter nebular gas us ing the theoretical Bry/BrQ- flux ratio of 0.35 
(IOsterbrocklll989l) . For M4, M5 and M6, we obtain NIR extinc- 
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tions of Av=13.5, 3.2 and 8.5, respectively. We use these values 
to de-redden the NIR line fluxes. 



4.2. Age and mass 

Together with Bry and Bra, we can use the measurements of 
the MIR fine structure lines [ArIII]8.9yum, [SIV] 10.4 yum and 
[Nell] 12. 8 /im in order to derive the cluster parameters. 
Predicting the intensities of the fine structure lines is more com- 
plex than in the case of the hydrogen lines. Their ionisation 
potentials are larger (respectively 27.6eV, 34.8eV and 21.6eV), 
which implies that not only the H ionising luminosity, QH, must 
be taken into account, but also that the shape of the ionising con- 
tinuum must be known or assumed. Moreover, the fluxes in these 
lines depend on the ionisation factor, closely linked to the ge- 
ometry of the source while we do not have much handle on it. 
Finally, the line fluxes depend on the element abundance ratios: 
given the location of the YMCs in their host galaxy, a set of solar 
abundances looks adequate even though a slight over abundance 
(up to a factor two) could be expected. 

What is remarkable on the spectra of the three sources, M4, 
M5 and M6, is that while the [Nell] line is conspicuous, only 
upper limits can be measured for the [Arlll] and [SIV] lines: 
the [NeII]/[SIV] line ratios are greater than 20, 100 and 140 for 
M4, M5 and M6 respectively, and the [NeII]/[ArIII] line ratios 
are larger than 10, 50 and 60. Such ratios are impossible to re- 
produce with very young stellar populations, for which the line 
ratio [NeII]/[SIV] is usually found to be lower than unity (see 
the SWS spec trum of R136 and for instanc e the spectrum of 
NGC5253 in iMartm-Hernandez et ani2005h . So, a substantial 
difference in the line ratios is unveiled, by one to two orders of 
magnitudes. 

In order to examine this differ ence, we use the CLOUDY 
emission line libraries computed b y Panuzzo et al] (20030 for 
branching into GRASIL dSilva et al. l ll998l) . In these libraries, the 
ionisation factor is varied by changing the gas density and the 
gas filling factor. iPanuzzo et alj (l2003l) demonstrate that differ- 
ent geometries with a similar final ionisation factor will produce 
very similar line fluxes. We compute the Bry, Bra, [Arlll], [SIV] 
and [Nell] line fluxes for a "cluster model" at ages between 1 and 
lOMyr with an "instantaneous" star formation history. As men- 
tioned above, we adopt solar abundances and we notice that in 
no way an overabundance by a factor two could account for the 
very large observed differences in line ratios. 

Since the [ArIII]/[NeII] line ratios observed in M4, M5 and 
M6 are very small, we use the libraries predicting the small- 
est values for this ratio: filling factor of 10"^ and density of 
10'*cm"''. The evolution with age of the un-extincted line fluxes, 
at the distance of NGC1365, are displayed in Fig. [15] for a total 
mass of stars of lO^M©. At first order, these line fluxes show, 
for reasonable values, little dependence on the amount of gas, in 
agreement with the ionisation-bound situation. 

In Fig. [161 we present a comparison between the modelled 
line fluxes and the observed ones. The procedure for building 
this figure for the three YMCs in NGC1365 was the following: 

- the extincted modelled line fluxes are computed, for each 
YMC, with the related extinction value derived above, 

- the extincted modelled line fluxes are then normalised, by 
changing the cluster mass, so that the modelled Bry line flux 
equals the observed Bry flux, 

- the line ratios between observed and modelled line fluxes 
are then computed for each age and plotted (left Y-axis) as a 
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function of age. Whenever the modelled flux falls inside the 
observational error-bar or is lower than the detection upper 
limit, then the ratio is set to unity, 
- for each age, the corresponding stellar mass is also plotted 
(right Y-axis). 

The model can be considered to "fit" the data if all three 
ratios are equal to one. This is absolutely not possible before 
6 Myr, where we would expect intense [Arlll] and [SIV] lines, 
one order of magnitude brighter than allowed by the upper lim- 
its measured on the spectra of the three YMCs. A reasonable 
"fit" can be obtained only for ages greater than 6 Myr. Notice 
that in the case of M5, there remains a small discrepancy, by a 
factor two, for the [Nell] flux, which, given the large differences 
involved and the first-order matches performed, is left aside. 

In Fig. [17] we plot the Starburst99 spectra of an instanta- 
neously formed Salpeter star cluster at ages 4, 5, 6 and 7 Myr, 
and compare them with the ionisation energy of [Nell], [Arlll] 
and [SIV]. This figure shows , between 5 and 6 Myr, a huge drop 
of luminosity (of about two orders of magnitude) for the con- 
tinuum in the region of the ionisation energies for [Arlll] and 
[SIV]. On the contrary, the continuum around the ionisation en- 
ergy for [Nell] only suffers a modest decrease. This shows that 
the virtual absence of [Arlll] and [Nell] lines after 6 Myr is due 
to the "absence" of hard enough continuum for these ion species. 

The firm conclusion to be drawn from Fig. [16] is a strong 
evidence for the YMCs in NGC1365 to be older than 6 Myr 
Because of the important fading suffered by stellar clusters along 
the first Myrs of their evolution, the fact that they are relatively 
old also implies that they are very massive. Let us consider the 
age of 7 Myrs: the computed stellar masses in this case for M4, 
M5 and M6 are then respectively 1.6,1.5 and 1.0 xlO^ M© for a 
IMq lower mass limit of the IMF. 

We show in Fig. [18] the stellar mass function in M4 with 
these parameters (age of 7 Myrs and mass of 1.6xlO^M0). The 
most massive stars in the clusters, 20-25Mo, are a few several 
lO'' in number, while low mass stars are expected to be present 
by millions. 

Consideration of a line flux library with lower filling factors 
shows that, while the [Arlll] line flux is not affected, the [SIV] 
line flux may increase by a factor up to 10. This remains com- 
patible with the data as long as the clusters are only about 1 Myr 
older than the age derived for a filling factor of 10"^. Hence, a 
safe range for the age of the YMCs is 6-8 Myr. We exclude the 
clusters to be older than 8 Myr since if it was the case, the masses 
derived would be of the order of IO'^Mq for the only stellar com- 
ponent and become inc ompatible with the CO measurements by 
ISakamotoet al.l(l2007h . 

This is consistent with the detection of the CO absorption 
lines at 2.3 fim in the spectra of M4, M5 and M6, sinc e this CO 
feature only appears after 6 Myr (iLeitherer et al.|[T999h . 

In principle, one could also use the ratio Hel 2.06 /jm/Bry to 
derive information about the hardness of the ionising continuum 
radiation, hence the age of the stell ar popula tion (DoyonetaTj 
119921) . However,(ShieMs| (Il993h and lEumsden et al.1 (l2003h have 
demonstrated that the use of this ratio is not trustworthy. 



4.3. Ionising photon emission rate 

Using the NIR extinctions derived above, we can compute the 
de-reddened hydrogen line fluxes, and get a direct handle on the 
ionising photon emission rates of the YMCs. For an ionisation- 
bound situation, justified for an embedded cluster, all the ionis- 
ing photons (A < 912A) emitted by the stars are absorbed. The 
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number of H ionising photons is then directl y proportional to th e 
flux in any specific recombination line. From lOsterbrodd (1 1 989h : 

Q[H+] = ^ X ^ , where ffB/<, ~ 2.96 (1) 

The Bry de-reddened fluxes of M4, M5 and M6 being re- 
spectively in the ranges [8.2-12.3], [7.6-11.5] and [13.7-20.6] 
xlO"'^ erg s"' cm"^, this leads to ionising photon emission rates 
in the ranges [3.4-5.1], [3.2-4.8] and [5.7-8.6] xlO^V. 

Let us check now whether the ionising photon emission 
rates derived from the Bry de-reddened fluxes of M4, M5 and 
M6 are consistent with the observed radio data in the cen- 
timetre range, given the ages just derived for M4, M5 and 
M6. The measuremen ts reported in Table [1] and performed by 
ISandqvist et alj d 19951) indicate negative indexes revealing, as 
expected for young clusters of the derived ages, the presence 
of a non-thermal component. The measurements also show that 
the indexes become increasingly positive with increasing wave- 
length (ffgj?™ > ffacm)- signature of optically thick 
radio emission. Indeed, it is known that UDHll regions can be 
optically thick even a t 6cm (|&)bulnickv & Johnson 199^. In 
ISandqvist et alj (Il995h . even though the 2 cm fluxes are not given 
explicitly, one can extrapolate their values using the 6cm flux 
and the 6cm/2cm spectral indexes (given in their Table 2). The 
2cm flux is interesting, since it is both optically thin and likely 
to be free from a non-thermal contribution. Under this assump- 
tion, the measured radio fluxes provide independent estimates of 
the io nising photon emission rates of the clusters. From Table 
2 from lSandqvist et al] d 19951) . we get the following 2 cm fluxes 
for M4, M5 and M6: 1.75 mJy, 0.61mJy, 1.70mJy respectively. 
The corresponding figures for the ionising photon emission rates 
are 2.0xl0^2s-', 7.0x10^' s^'and 1.9xl0^2s-'. These values 
are slightly lower than, but in reasonable agreement with those 
derived from the de-reddened emission lines. 

4.4. Final remarks about the derived YMC parameters 

Along this procedure, the derivation of the cluster mass and the 
derivation of the ionising photon emission rate do depend on the 
validity of the assumption of a foreground uniform extinction. 
Conversely, the age derivation is quite robust and relies essen- 
tially on the disappearance of some emission lines relatively to 
others. The factor to be matched is huge, one to two orders of 
magnitude, and can in no way be understood in terms of mea- 
surement uncertainties. 

Would there be a possibility to explain the observed line 
ratios, maintaining a very small age (around 1 Myr)? Indeed, 
instead of the Salpeter IMF which we consider, we might in- 
voke an ad-hoc fancy IMF with a low upper cutoff (at about 
23Mo), but such a value is not justified. In addition, the ob- 
served negative spectral indexes of the radio emission from the 
YMCs constitute an independent argument favouring advanced 
ages, since the non-thermal component of the centimetre emis- 
sion traces the presence of supernovae, which only occur after 
3-4 Myr. Therefore, we retain our age determinations as provid- 
ing the most sensible interpretation. 

For the quoted age range for the clusters (6-8 Myr), the cor- 
responding stellar masses for the three clusters are in the range 
[0.6-3.]xlO^Mo. These masses are derived from the unextincted 
Bry fluxes and a Salpeter IMF with a lower boundary at 1 M©. 
This masses must be multiplied by 2.5 if the lower boundary of 
the IMF is decreased to 0. 1 Myr. 

The uncertainty on the estimate of the deredenned line lumi- 
nosities , which we recall is performed under the assumption of 



foreground extinction, is of the order of only 20%, hence small 
with respect to the uncertainty due to the uncertainties on the age 
and the IMF shape. 

We hence conclude that the clusters are 6-8 Myr old, with 
masses of the order 10^ Mq and ionising photon emission rates 
of several 10^^ s"'. 

All these elements concur to the stunning result that the three 
MIR/radio sources in NGC1365 are among the most massive 
clusters observed so far. Surprisingly for their relatively "ad- 
vanced" age, their MIR emission indicates that they have not 
swept away all of their surrounding material yet. Usually, clus- 
ters are believed to remove very quickly the material in which 
they formed, within a few million years. The peculiar case of 
these YMCs opens the question of whether their extreme mass 
and their location close to the dust lane of their host galaxy play 
a role in the fact that, at 7Myrs, they still contain important 
amounts of gas. 

5. Fitting the YIVIC NIR/IVIIR SEDs 
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Fig. 17. Starburst99 model spectra for an instantaneously formed 
Salpeter cluster at ages 4, 5, 6 and 7 Myr The vertical lines show 
the first ionisation potentials forH, [Nell], He, [Arlll] and [SIV]. 



In this section, we discuss the fitting of the NIR/MIR SED of 
the YMCs. We later add the constraints from thek emission line 
spectrum. 

5.1. Analysis of the NIR/MIR SEDs through comparison with 
R136 

First, we examine in Fig. [19] a comparison of the SEDs of M4, 
M5 and M6 with the ISO CVF spectrum of R136, scaled to 
the distance of NGC1365, normalised at 12.5 yum, and subject 
to some foreground extinction: we use a set of foreground ex- 
tinctions in the range Av=0 to Av=40. 

One sees that in the 8-9 //m window, the three NGC1365 
YMCs are brighter than R136 by factors between 2 and 10. This 
factor cannot be just explained by a calibration error for the 8- 
9 fim range (See Sec. I2.4l i since, in the worst case, we evaluated 
a flux uncertainty by a factor 2 for this region of the spectrum. 
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Fig. 20. Comparison of the observed spectra of M4, M5 and M6 (black) to the scaled ISO CVF spectrum of R 136 (brown) and to 
the predicted spectra of the 2-component model described in the text. For the modelled spectra, the colour coding is the following: 
purple is the total spectrum (properly scaled for comparison with the source spectra); red is the total spectrum of the molecular cloud 
component (the optically thick one) component; blue is the spectrum of the dust (essentially very small grains, VSG) emission of 
the ciiTus component (the optically thin one) and green is the spectrum of the[stellar+gas] component of the ciiTus component 



Even though a full N-band spectrum would be mandatory to de- 
cisively settle this difference, let us put forward the possibility 
that the high 8-9 /vm continuum is due to the presence in the 
YMCs of a bright 8.6 /im PAH feature. In the usual PAH tem- 
plates (as used in GRASIL, iSilva et al.]| 1998b . the 12.7 fim fea- 
ture is much brighter than the 8.6 /zm feature. However, the 
8.6)Um feature can be significantly boosted for a high degree 
of ionisation of the PAH molecules. Indeed, for ionised PAHs, 
the ratio of cross sections cr%,^lcr\2j is enhanced by a factor 7-8 



with respect to the neutral case jLi & Drainell2001l;lDraine & Lil 
l2007h . 

About the 10-13/im part of the SED, its slope can be rea- 
sonably matched by a foreground extincted version of the R136 
spectrum. In particular the lO/im-13 //m VISIR spectra are well 
reproduced for M5 and M6 (including the 11.3 yum PAH feature). 
For M4, the resemblance is not so obvious, but the M4 SED is 
also noisier. 
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Fig. 18. Mass function for M4, given the age and mass derived 
in Sec. 14:21 




Av= 0, 10, 20, 30, 40 



_i I 1 1 I 1 1 1 I 1 1 1 1_ 



wavelength [^m] 



10 12 14 




wavelength [/im] 



Fig. 19. Comparison between the observed spectra of 
NGC1365/M4, M5 and M6 to the ISO CVF spectra of 
R136 suffering from distinct foreground extinctions, from Av=0 
to Av=40. The CVF spectra are normaUsed to the VISIR spectra 
at 12.5 ;um. 



M4: age mc=4.8Myr, Mcirrus = 8.4e + 06Msol, Mmc=1.9e + 06Msol 
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Fig. 21. Quality of the 9 predicted observables for the 2- 
component model for M4, M5 and M6 (see Sec. |5j. The ob- 
served uncertainty is shown with the grey stripe. If obs„,i„ and 
obs„,ax represent the lower and upper limits of the error bar, then 
the value on the Y-axis is {mod-obs,„i„)/(obs,„ax-obs,„i„), where 
mod is the predicted value. For the observables for which we 
only measure an upper limit, then a value Y=0 is given to the 
corresponding point if the mode led value is lower than this limit. 
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Fig. 22. Threshold mechanical lumi nosity for two heatin g ef- 
ficiency (reproduction of Fig. 4,b of IWiinsch et alj|2007[) . The 
horizontal line shows that, according to our estimates, the three 
MIR/radio clusters of NGC1365 lie in the bimodal regime re- 
gion. 
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In summary, the SED of M5 resembles that of R136, ex- 
cept in the 8 //m-9 fim region where ionised PAH, boosting the 
8.6jum feature, could be an explanation. Therefore, M5 could 
be seen as a scaled-up version of R136, from the SED point of 
view. However, the emission line ratios discussed earlier clearly 
show that M5 cannot have the same age as R136, but is at a 
more advanced stage 6-8 Myrs. This is not a contradiction be- 
cause the SED matching is not an age estimator The key issue 
is to understand how M5, at 7 Myrs, could have retained some 
dust component? 

What is the conclusion of the comparison with R136 for M4 
and M6? A very distinctive feature observed in M4 and M6 is 
their steep continuum rise between 3 jum and 4 jum , in contrary to 
the spectrum of R136. This type of feature is reminiscent of the 
emission from an optically thick dust component and suggests 
its presence in the YMCs. Therefore, we investigate hereafter 
such a possibility, using a more complex modelling tool. 

5.2. Modelling of the NIR/MIR SEDs 

Let us compare the SEDs we have in hand for the three YMCs 
to the SED prediction of a cluster including two components: (a) 
an optically thick component, which we will call the "molecu- 
lar cloud" component, with stellar mass M„,c [which we would 
assume to be associated with the "youngest" stellar population 
in the cluster, in case the star formation would not have been in- 
stantaneous] and (b) an optically thin component, which we will 
call the "cirrus" component, with a stellar mass Mdrms [which 
we assume to be associated with the gas and with the "oldest" 

stellar population in the cluster] . 

We use GRASIL dSilva et al.lll998h to model the SED clus- 
ter. In a broad framework, GRASIL is designed to model the 
spectro-photometric evolution of dusty galaxies, including a de- 
tailed treatment of the radiative transfer through dust. It can be 
used to simulate the evolution of the properties and emission 
of a stellar population with the two components as described 
above: (a) the "molecular cloud component" for which the stars 
are still located at the centre of their parent molecular cloud and 
for which a large optical thickness is expected, and (b) the "cir- 
rus component" in which both the stars and the dusty gas show 
an extended distribution, and for which the optical thickness is 
globally low. 

Fig. |20] illustrates how such a two-component model can 
match the global features of the observed SEDs of the three 
YMCs. This figure displays the observed SEDs for M4, M5 and 
M6 (in black). For each source, the R136 spectrum scaled at 
around 2 fim is superimposed (in brown). As well, the different 
components of the model are superimposed and shown under the 
following colour code: in green is the joint stellar emission and 
continuum emission from the HII region arising from the "cir- 
rus" component, in blue is the dust emission from the "cirrus" 
component, in red is the net emission from the thick "molecu- 
lar cloud" component (dust and stellar emission, but essentially 
dominated by the dust) and in purple the total emission. The key 
parameter for the modelled "molecular cloud" component is its 
optical depth. For each YMC, this optical depth has been cho- 
sen in order to fit the slope in the red wing of the silicate feature 
around 10 /im, observed in the YMC, and has been derived after 
full treatment of the radiation transfer. The extinction values Av 
which are required to reproduce the observed silicate feature red 
wing slopes are very large: respectively 100, 70 and 130, for M4, 
M5 and M6. Indeed, for the modelled "molecular cloud" com- 
ponent, dust is simultaneously the source of emission and ex- 
tinction. For details about the radiative transfer in GRASIL, see 



ISilva et aTl ( Il998l) . Notice that the eff'ect of such an extinction 
cannot be compared in a straightforward manner to the effect of 
the simple foreground extinction shown in Fig. [19] at the begin- 
ning of this section, and estimated by comparing the YMC SED 
to the R136 spectrum. 

Following this modelling, the NIR 3-4 fim slope and the MIR 
slope between lOyum and 13 //m are both steep as a result of the 
emission of the "molecular cloud" component and an important 
bump on the SED in the 4-8 jum range is predicted. This is con- 
sistent with the observed rise of the SED in the 3-4 fj.m range, in 
M4 and M6. 

In conclusion, the NIR/MIR modelled SEDs in Fig.l20lpro- 
vide qualitative matches for the observed SEDs in M4 and M6. 
In the case of M5, the slope of the continuum in the NIR 3- 
4 fim range is overestimated in the model and the presence of a 
"molecular cloud" component does not seem to be mandatory. 

Going one step further, let us investigate the possibility that 
the "molecular cloud" component represents a second stellar 
generation (hence younger), still deeply embedded in dust. For 
this purpose, we perform fits of the two-component model pre- 
dictions using the following observational constraints: the Bry, 
Bra, [Arlll], [SIV] and [Nell] fluxes, and the R, Ks, L' and 
12.5 fim flux densities. The age of the cirrus component is set 
to 7 Myr (derived earlier from the analysis of the line ratios) 
and the optical thickness of the "molecular cloud" component 
to the values given above and derived from the silicate red wing 
slope. The free parameters in the fits are the age and mass M„,c 
of the "molecular cloud" component (the age is forced to be 
< 6 Myr), the mass M„rn« of the "cirrus" component and the 
foreground extinction. Given the observational uncertainties and 
because the centre of the error-bar does not necessarily represent 
the most probable value of the observable, different solutions can 
be found. Yet, the results are globally identical. The quality of 
the fit for one solution is shown graphically in Fig. [21] and in 
fact. Fig. l20l show this same solution. For each observable, the 
error-bar has been scaled to the range [0-1]: the values of the 
modelled observables are given on this same scale. A dot inside 
the grey zone means that the modelled observable falls inside the 
error-bar 

The solutions shown in Fig. |2T| are good for M4 and M6, 
and marginally acceptable for M5. Indeed, Fig.|20]was akeady 
showing that the two-component model predicts a steep slope in 
the 3-4 fim range which is observed for M4 and M6, but not for 
M5. For the fits presented in Fig. [21] the age of the "molecu- 
lar component", the stellar mass of the "cirrus" component and 
the stellar mass of the "molecular cloud" component are indi- 
cated above each plot. For M4, the stellar mass in the "molecular 
component" is four time lower than in the "ciiTus" component. 
On the contrary, it is almost two times greater in the case of 
M6. Indeed, in Fig.|20]the relative importance of the "molecular 
component" (red) to the "cirrus component" (green and blue) is 
clearly greater in M6 than in M4. 

The two-component model presented here, even though sim- 
plifying a complex real situation and at the same time involving 
a large number of parameters, allows an interesting qualitative 
match to the observables. This is true to a good extent for M4 
and M6. Notice that the two-component model predicts a bump 
of the SED in the 4-8 jum range, differing from the R136-type 
SED. The SED of M5, on the contrary, looks more like that of 
R136, and does not require the presence of a thick "molecular 
cloud" component. This might be related to the fact that M5 
is seen lying slightly outside the dust lane, while M4 and M6 
project right in the dust lane. Indeed, the fact that the pressure of 
the surrounding interstellar medium is certainly higher around 
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M4 and M6 than around M5, may have affected their evolution 
and hence their intrinsic properties, and in peculiar their ability 
to remove the material from which they formed initially. 

5.3. Conclusions on the YMCs parameters 

Let us summarise the conclusions reached after the successive 
steps of the data analysis. 

First, we have found that the three YMCs in NGC1365 are 
compact clusters. This is a solid conclusion, coming from a di- 
rect comparison with the cluster R136 in the LMC. 

Second, the age estimates of the YMCs come from three in- 
dependent channels: 

(a) their NIR emission line ratios, with intense [Nell] and 
undetectable [Arlll] and [SIV] lines, can be understood only if 
their age is in the range 6-8 Myrs, 

(b) their radio indexes are steep, suggesting that the non- 
thermal contribution (supernovae remnants) is larger than the 
thermal contribution (HII region): although it is not possible to 
extract a precise value of the cluster age from its radio index, we 
know that the clusters must be older than 3 Myrs, 

(c) the detection of the CO absorption lines at 2.3 yum in their 
spectra tells us that the YMCs are older than 6 Myrs. 

Therefore, we find the YMC ages to be around 7 Myrs. This 
again is a fairly solid conclusion. 

Third, given their ages and their ionising photon emission 
rates (hence stellar luminosities), the stellar masses of the YMCs 
are found to be of the order of lO^M©, implying a mass of the 
initial molecular material from which they have been formed of 
at least several IO^Mq. The mo l ecular mass deduced from CO 
observations bv lSakamoto et al] (l2007h . of the order of lO'^M©, 
is indeed consistent with the figure we have obtained. 

Fourth, we find that, after 6 to 8 Myr of evolution, still impor- 
tant amounts of gaseous/dusty material, traced by the radio/MIR 
emission and nebular lines, are found within or in the vicinity of 
the star cluster. Since star clusters usually sweep away their gas 
on a much shorter time scale, our observations present a case for 
trapping of material in star clusters, at very high masses and in 
a dense environment. This gas trapping can lead to peculiar evo- 
lution of the star clusters by allowing subsequent star formation 
events to happen. 

In the following, we discuss, in the framework of the theo- 
retical model by Tenorio-Tagle and co-workers, how such a gas 
trapping can be explained. 



6. Gas trapping in the extremely massive clusters in 
NGC1365 

Theoretical piec es of work by the group of Silich and 
Tenorio-T agle (T enorio-Tagle et al.l |2007[ IWtinsch et all 120071: 
ISilich et al.l 120071) analyse in details the hydrodynamics of the 
matter re-inserted in young clusters by stellar winds and super- 
novae ejecta. They show that in the case of massive and compact 
clusters, the re-inserted material is exposed to strong radiative 
cooling in the central parts of the clusters, because of the high 
density of the gas in these regions. This strongly aff'ects the dy- 
namics of the gas and can lead to a bi-modal hydrodynamic so- 
lution: (i)) the matter injected inside a certain radius called the 
stagnation radius is accumulated and eventually becomes gravi- 
tationally unstable leading to further star formation and, (ii)) out- 
side the stagnation radius, the re-inserted material flows out of 
the cluster, building up a stationary wind dWiinsch et al.ll20()7l) . 



One candidate for such a bi-m o dal ga s behaviour has been 
identified in M82 by ' Sihch et al.1 ( l2007h . using observations 
from Smith et al. (2006). The observed parameters of the cluster, 
called M82-A1, imply such a bi-modal hydrodynamic solution. 
The adopted model leads to a much reduced rate of mass deposi- 
tion in the interstellar matter, and a much reduced wind terminal 
velocity, compared to the adiabatic wind model. 

The three clusters NGC1365/M4, M5 and M6 are also ex- 
cellent candidates for displaying such a bi-modal behaviour: 
they are both very massive and compact, hence the gas density 
in their central regions is high. Also, the fact that with an age 
7 Myr, these clusters have still retained an important amount of 
gas and dust, as shown by their MIR emission, suggests that the 
gas removal mechanism has not been eflicient. Finally, the sus- 
pected "molecular cloud" component which we derive from the 
SED modelling of M4 and M6 - the two YMCs which project 
onto the galaxy dust lane - could trace on-going star forming 
events, which recalls a prediction from the bi-modal model. 
Nevertheless, in the case of M6, we find that M^c > Mdnus- 
This is incompatible with the idea that, for M6, the "molecular 
component" traces an ongoing star formation involving only the 
material re-inserted by the "first" generation stars. The higher 
pressure of the dust lane gas in which these two clusters seem 
to be embedded could also have increased the efficiency of the 
recycling of the matter into new generations of stars. 

In order to investigate quantitatively whether the MIR/radio 
emitting clusters in NGC1365 are viable candidat es for this spe- 
cial re gime, we have used the results published in lWiinsch et al.l 
(|2007|) . Simple analytic formulae are provided in order to test if a 
cluster is undergoing such a bi-modal hydr odynamic solu t ion, or 
if it only drives a stationary wind. Fig. 4 in lWiinsch et al.l(l2007h 
displays the threshold mechanical luminosity as a function of the 
star cluster radius, that separates single stationary wind solution 
from bi-modal solutions. The threshold is shown to depend on 
the heating efficiency (parameter 77 in their work). This parame- 
ter represents the efficiency of the thermalization of the gas. 

To estimate the mechan ical luminosity of th e clusters, we 
use the Starburst99 models (Leitherer et al.|[T999h : for a Salpeter 
IO^Mq model, with solar metallicity, for an age range be- 
tween 5 and 9 Myr, the range of mechanical luminosities is 
lO'"'^ to lO"^^'^ erg s"'. The radii of the cluster stellar component 
are lower than lOpc. This radius is directly measured for M5, 
and assumed for M4 and M6, where no clear visi ble counterpart 
is dete cted. In the radius vs. luminosity diagram of lWtinsch et al.l 
(|2007|) . reproduced in Fig. |22] the clusters clearly fall in the bi- 
modal solution region. 

The picture drawn here of the three bright MIR/radio sources 
seems coherent. In the massive reservoir of matter surrounding 
the nucleus of NGC1365, and shown by the prominent dust lane, 
very massive and compact clusters have formed. The hydrody- 
namics of the gas inside these clusters (re-inserted and/or unused 
for star formation) is susceptible to follow the bi-modal hydro- 
dynamic solution described above. This leads to a much less effi- 
cient removal of the gas, which stagnates in the inner regions of 
the clusters, while it is ejected in the outer regions with reduced 
velocities, compared to common adiabatic cluster outflows. This 
could be traced by the velocity gradient shown in Fig. |9] but we 
leave this analysis for a future publication. 



7. Conclusion 

Let us summarise the main conclusions we have reached along 
this work. 
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The three MIR/radio sources in the circumnuclear star form- 
ing ring in NGC1365 have been found to be very compact and 
very bright embedded clusters. Under the simple assumption of a 
foreground extinction, their MIR emission line spectrum points 
towards ages of about 7 Myrs. Such an age is also consistent with 
other age indicators like the slope of their radio centimetre emis- 
sion or the occurrence of thir NIR spectrum of the CO absorp- 
tion feature. Given their ionising photon emission rate and age, 
these YMCs are found to be extremely massive objects, with a 
mass of at least 10^ M©. They are the most massive star clus- 
ters ever found so far. They must have formed from molecular 
clouds of several 10^ M©. This figure is consistent with the pres- 
ence within the NGC1365 central 2kpc diameter region, where 
they are found, of large amounts of molecular gas, of the order 
of IO^Mq, as derived from the CO molecule emission. 

The remarkable output of our analysis of the data is the fact 
that such embedded YMCs have retained an important amount 
of gas and dust, in spite of their relatively advanced age, around 
7 Myrs. This a quite puzzling and interesting fact, which may be 
related to the extreme mass of these three YMCS. Indeed, some 
theoretical works have analysed the fate of re-injected gas in 
very massive clusters and found that the gas may be trapped and 
even lead to secondary star forming events. The three MIR/radio 
clusters in NGC1365 are good candidates to test the predictions 
of such models. 
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